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We perform a global fit of the parameters of the Standard Model with a sequential fourth gen- 
eration (SM4) to LHC and Tevatron Higgs data and electroweak precision data. Using several 
likelihood ratio tests we compare the performance of the SM4 and SM3 at describing the measured 
data. Since the SMS and SM4 are not nested (i.e. the SMS can not be considered as a special case of 
the SM4 with some parameters fixed) the usual analytical formulae for p-values in likelihood ratio 
tests do not hold. We thus apply a new method to compute these p-values. For a Higgs mass of 
126.5 GeV and fourth-generation quark masses above 600 GeV we find that the SM4 is excluded at 
S.lcr. 



METHOD AND INPUTS 

In this paper we study the SM4, which differs from the 
established Standard Model (denoted by SMS) by an ad- 
ditional fermion generation. We treat the masses of the 
extra fermions as free parameters and allow for arbitrary 
flavor mixings among the quarks of the four generations 
in our fits. Large mixings of the fourth-generation lep- 
ton doublet with those of the first three generations are 
ruled out [1 from data on lepton-fiavor violating decays 
and lepton-fiavor universality [2 . Recent NA62 data con- 
strain these mixing angles even further [3] . Including lep- 
ton mixing within the allowed range has a negligible im- 
pact on the electroweak precision observables (EWPOs). 
In the absence of lepton mixing the decay of the Higgs 
boson into neutrinos is invisible as long as the fourth- 
generation charged lepton is heavier than the correspond- 
ing neutrino. This invisible Higgs decay mode increases 
the total Higgs width and potentially counterbalances the 
effect of the enhanced gg ^ H production mechanism 
[H IS] , because the branching fractions into the observed 
final states are reduced [6l . Allowing for (even small) 
mixing of the fourth with the other lepton doublets can 
render the neutrino decay mode visible. Since we want 
to quantify the level at which the SM4 is ruled out, we 
may confine ourselves to the most conservative scenario 
with an unmixed fourth-generation lepton doublet. Like 
the SM3, the SM4 can be studied with Dirac or Majo- 
rana neutrinos. In the fits presented in this paper we use 
Dirac neutrinos. In our conclusions we briefly discuss 
the (marginal) changes in the results expected for Ma- 
jorana neutrinos. From a model-building point of view, 
the hierarchy between three almost massless neutrinos 
and a fourth neutrino with mass of order of the elec- 



troweak scale can be motivated by a symmetry enforc- 
ing massless neutrinos in the exact symmetry limit: e.g. 
three right-handed neutrino fields might carry some U(l) 
charge while the fourth neutrino field and the left-handed 
lepton doublets are uncharged under this new symmetry. 
The Yukawa couplings are small spurious breaking this 
symmetry, leading to three tiny neutrino masses and tiny 
mixings between the fourth and the other generations. 

A sequential fourth generation of fermions decouples 
neither from the production cross section cr{gg — )■ H) nor 
from the Higgs decay rate into photons. Consequently, 
current LHC Higgs data put the SM4 under serious pres- 
sure [7l410j. In a recent publication [7' we presented a 
global fit of the SM4 parameters to EWPOs and Higgs 
signal strengths measured at Tevatron and the LHC. The 
signal strength is defined as 

M^-^-n = 4^^^41f^^#^ . (1) 
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Here we update our results with all available data and 
analyse the status of the SM4 prior to the ICHEP2012 
conference. We also compute the statistical significance 
(p-value) at which the SM4 is excluded. As explained 
in U the computation of the p-value is non-trivial: due 
to the non-decoupling nature of the fourth-generation 
fermions the SMS can not be regarded as a special case 
of the SM4, i.e. the two models are not nested. Analyti- 
cal formulae for p- values only hold for nested models and 
thus the p- value of the SM4 has to be computed numeri- 
cally. To this end, a new CH — h framework for maximum 
likelihood fits and likelihood ratio tests called mi/Fitter 
was written. The implementation is discussed in |12j . 
In total, the following aspects of our previous analysis 
have been improved: 
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1. The masses of all four fourth-generation fermions 
are now consistently treated as free parameters. To 
avoid non-perturbative Yukawa couplings and con- 
straints from direct searches of fourth-generation 
quarks we require 600 GeV < mf , nib' < 800 GeV. 
We are aware that for fermion masses of 800 GeV 
the validity of perturbation theory is questionable 
at best. However, reducing the upper limit for the 
fermion masses can only lead to larger values 
in the SM4 and thus to smaller p-values. In this 
sense, the upper limit of 800 GeV is a conservative 
estimate. 

2. The signal strength for pp ^ H ^ tt measured at 
the LHC [13] is included in the analysis. 

3. In the global fit, the Higgs mass is no longer fixed at 
125 GeV, but is allowed to float in the range where 
experimental data on the Higgs signal strengths is 
available, i.e. 115 GeV < tuh < 150GeV.[Tl] 

4. Since, for a variable Higgs mass, no separate H 
77 signal strengths for the gluon fusion and vec- 
tor boson fusion production modes are available 
we only use the combined signal strength for pp — > 

— >■ 77 as input. 

5. For the two cases tuh = 126.5 GeV (the preferred 
Higgs mass of the SM3) and mn = 147 GeV (the 
preferred Higgs mass of the SM4) we perform likeli- 
hood ratio tests to compare the performance of the 
SM3 and SM4 at describing the measured data. 

Regarding the last point, a few more comments are in 
order. In likelihood ratio tests the difference Ax^ of min- 
imal values obtained in the SM3 and the SM4 is used 
as a test statistic. One then assumes that the measured 
observables are random variables distributed around the 
prediction of one model (e.g. the SM4) with a spread de- 
termined by their errors and computes the probability 
(p-value) that a random set of "toy-observables" leads 
to a Ax^ which is more extreme (e.g. more SM3-like) 
than the Ax^-value obtained from the real data. Note 
that this is different from the goodness-of-fit analysis pre- 
sented in [Hj, which used the x^ value of the SM4 as a 
test statistic and therefore did not compare the perfor- 
mance of the SM3 and the SM4. Also, the H ^ tt signal 
strengths were not included in their analysis. 

Unfortunately, the likelihood ratio tests can not be 
done (by us) if the Higgs mass is treated as a free pa- 
rameter. In that case, the signal strengths measured in 
each invariant mass bin of each Higgs decay mode would 
have to be treated as separate observables, and we do not 
have any information on statistical correlations between 
adjacent bins. Thus we only perform likelihood ratio 
tests for specialisations of the SM3 and SM4, where the 
Higgs mass is fixed to mH = 126.5 GeV (the value pre- 
ferred by the global SM3 fit) or niH = 147 GeV (the value 



preferred by the global SM4 fit). Then only the signal 
strengths at mn = 126.5 GeV and ruH = 147 GeV have 
to be treated as independent observables and correlations 
between these observables can safely be neglected. 

Note, however, that the information from all invariant 
mass bins is encoded in our x^ function. So, for example, 
the x^ value at mn — 147 GeV has a contribution due 
to the fact that there is a signal at itlh = 126.5 GeV. If 
the model under consideration had a Higgs boson outside 
the discovery reach of LHC (or no Higgs boson at all), 
the theory prediction for all signal strengths in all invari- 
ant mass bins would be zero. This leads to a constant 
contribution to the x^j which we are allowed to drop. 
Now assume that the model has a Higgs boson with some 
mass niH and a predicted signal strength /ith(™_ff)- Let 
fj-cximn) and ApL{mH) be the measured signal strength 
and experimental error for the corresponding invariant 
mass bin. After dropping the constant, the x^ function 
is 

2/ N [i^&{rnH) - j^cx{rnH)f -[jleximH)? 
X [mH) = ^r^TT Tf^ • (2) 

If there is a clear signal at the Higgs mass mjy , the second 
term gives a large negative contribution to the x^ func- 
tion. This contribution is not present if mn is in a region 
without a signal, so the minimum of the x^ function will 
usually be at a Higgs mass close to the signal. 

In the present analysis, the following experimental in- 
puts are used: 

i) fiipp -^H-^ WW*) measured by ATLAS [H], 

ii) fi{pp H ^ jj) measured by ATLAS [H], 

iii) fi{pp -> HV Vbb) measured by CDF and DO 

iv) fj.{pp ZZ*) and /t(pp — > _ff — > rf) mea- 
sured by ATLAS [11], 

iv) the electroweak precision observables (EWPOs) 

Mz, Tz, (Jhad, ^FB' ^FB' ^FB' "^c, ^6, Rl — 

Fj+z-ZPhad, Rc, Rb, sin^ measured at LEP and 
SLC [18 as well as toj, Mw, Fw and Aa{,^2i I^- 

v) the lower bounds mt'.b' 600 GeV(from the LHC) 
[20H23] and m/^ > 101 GeV (from LEP2) [19 . 

Unfortunately, there is no data for signal strengths as a 
function of the Higgs mass from CMS. 

On the theory side, the global fits with a variable Higgs 
mass were done with the CKMfitter software [23]. The 
EWPOs in the SM4 were calculated with the method de- 
scribed in [25], using FeynArts, FormCalc and LoopTools 
[26H28] to compute the SM4 corrections to the EWPOs. 
The EWPOs in the SM3 were calculated with the ZFit- 
ter software [^51 - I5T| . The Higgs width and branching ra- 
tios in the SM4 and SMS were calculated with HDECAY 
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FIG. 1. Minimum x values for a fixed neutrino mass as a 
function of m^^. The blue (red) lines show the results from 
the combined analysis of EWPOs and Higgs signal strengths 
with (without) the H ^ tt channel. The solid and dashed 
lines correspond to the SM4 and SMS, respectively. 



V. 4.45 [32], which implements results of |33H36] . The 
SM3 Higgs production cross sections were taken from 
[57] (LHC) and [511131] (Tevatron). For the numerical 
integration required to compute the p-values we use the 
Dvegas code [JD] which was developed in the context of 

[Hiia. 



RESULTS 

To show the impact of the 7? — >■ rr signal strength we 
plot the minimal value with and without the H — 
TT input as a function of the mass m^^ of the fourth- 
generation neutrino in Fig. [l] We see that for < 
60 GeV the minimum values are almost the same with 
and without the H ^ tt input. For rriu^ > 65 GeV the 
H ^ TT input increases the minimum by more than 
20. We also see that without the H ^ tt input the 
SM4 favours large values of .With the H ^ tt signal 
strengths included, the smallest values are obtained 
for m^^ between 50 and 60 GeV. 

This can be understood as follows: the production rate 
of Higgs bosons in gluon fusion is enhanced by a fac- 
tor of 9 in the SM4 due to the contributions from addi- 
tional heavy quark loops. On the other hand, the effec- 
tive HWW, HZZ and H^^ couplings are suppressed by 
the higher order corrections discussed in [36]. No such 
suppression is possible for H tt, so we would expect 
a — > TT signal strength of 9. The only way to reduce 
this signal strength is to open the invisible H — >■ 1/41/4 
decay mode, which then suppresses all branching ratios 
by a common factor. Thus, for large values of m,^^, the 
fit gets considerably worse if the H tt channel is in- 
cluded. 

Figs. [2] and [3] show the minimum value as a function 
of the Higgs mass mn in the SMS and SM4, respectively. 
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FIG. 2. The minimum value of the SMS as a function of 
the Higgs mass mn- The solid line shows the results of the 
combined analysis of signal strengths and EWPOs. For the 
dashed line only the signal strengths were included in the fit. 



The solid lines show the results of the combined analy- 
sis of signal strengths and EWPOs while for the dashed 
lines only the Higgs signal strengths (including H — tt) 
were used as inputs. We see that the SM3 clearly prefers 
a Higgs mass near 126.5 GeV. This is in agreement with 
a similar analysis presented in [43] . There is another 
local minimum at mu = 145 GeV, but with a consid- 
erably larger value. The x^ function of the SM4 in 
the combined analysis of signal strengths and EWPOs 
also has one minimum at ran = 126.5 GeV and another 
one at tuh = 147 GeV. Here, the x^ values are almost 
the same, but still larger than the minimal x^ value of 
the SM3, x^ = 23.3, by about 8 units. Note that for 
non-nested models or models with bounded parameters 
the relation between x^ values and p-values is no longer 
given by Wilks' theorem. Thus, in the case of the SM4, 
the number of degrees of freedom is an ill-defined con- 
cept and the p-values have to be calculated by numerical 
simulation. For the simulations we used the mj/Fitter 
package [TT]. Further details on the statistical issues 
and the myFitter simulation method can be found in 
[T^ . For niH = 147 GeV the signals at invariant masses 
near 126.5 GeV would be interpreted as statistical fluctu- 
ations. Then the data would be better described by the 
SM4 because it has more mechanisms for suppressing its 
Higgs signals. These mechanisms were discussed in [7]. 

Fig. |4| shows the pulls of the Higgs signal strengths 
for the SM3 with a Higgs mass of 126.5 GeV and the 
SM4 with a Higgs mass of 126.5 GeV or 147 GeV. We 
see that in the SM4 with ttlh — 126.5 GeV the measured 
H ^ TT signal strength deviates by more than 4(t from 
its predicted value. This is due to the effect mentioned 
in the discussion of Fig. [ij For the SM4 with ran — 
147 GeV the measured signal strengths for the invariant 
mass bin at 147 GeV are in good agreement with their 
theory predictions. However, in that case the x^ receives 
a large contribution due to the fact that the measured 
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TABLE I. p-values obtained from the likelihood ratio tests 
for fixed Ifiggs mass. Both, the SM4 with niH ~ 126.5 GeV 
and 7nH = 147 GeV are compared to the SM3 with a fixed 
Higgs mass of 126.5 GeV. In the first row, only the Higgs 
signal strengths were used as inputs. The second row contains 
the results of the combined analysis of signal strengths and 
EWPOs. The number of standard deviations corresponding 
to each p- value are shown in parentheses. 



FIG. 3. The minimum value of the SM4 as a function of 
the Higgs mass itih. The solid line shows the results of the 
combined analysis of signal strengths and EWPOs. For the 
dashed line only the signal strengths were included in the fit. 
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FIG. 4. Pulls of the Higgs signal strengths for the SMS with 
niH = 126.5 GeV and the SM4 with mn = 126.5 GeV and 
niH = 147 GeV. 



values of the signal strengths in the invariant mass bin 
at 126.5 GeV deviate from their predicted values of zero. 

Table |l] shows the p-values obtained from the likeli- 
hood ratio tests for the two SM4 Higgs masses. We see 
that, based on the Higgs signal strengths alone, the SM4 
scenario with mn — 126.5 GeV is ruled out at almost 
4(7 while the scenario with tuh = 147 GeV is only ex- 
cluded at 3a. At a fixed Higgs mass of 126.5 GeV the 
electroweak fit is actually better in the SM4 than in the 
SMS. Thus, if the EWPOs are included in the fit, the p- 
value increases to 2 permille, which corresponds to S.lcr. 
The lower bound mf^b' ?J 600 GeV is not essential for 
this result, relaxing this bound to mt',b' ^ 400 GeV de- 
creases the minimum-x^ by 0.6. For the SM4 scenario 
with iriH — 147 GeV the p-vahie drops to 0.74 permille 
(3.40-). In any case, the SM4 is excluded at more than 
3cr. 



CONCLUSIONS 

We presented a combined analysis of Higgs signal 
strengths and EWPOs in the context of the Standard 
Model with three or four fermion generations. The SM3 
is in good agreement with the experimental data and the 
best-fit Higgs mass is 126.5 GeV. The SM4, on the other 
hand, struggles to describe the Higgs signal strengths 
measured at Tevatron and the LHC. The function 
of the SM4 has two minima at niH = 126.5 GeV and 
niH = 147 GeV with essentially the same value, which 
is larger than the minimal x^ value of the SM3 by 8 units. 
The second minimum of the SM4 function occurs 
because the SM4 cannot reproduce the signal strengths 
measured at 126.5 GeV very well, so that an SM4 with a 
Higgs mass nowhere near the observed signals describes 
the data equally well as an SM4 with itih = 126.5 GeV. 
To quantitatively compare the performance of the SM3 
and SM4 at describing the data we performed likelihood 
ratio tests for fixed Higgs masses niH of 126.5 GeV in the 
SM3 and mH ^ 126.5 GeV, 147 GeV in the SM4. The 
p-values were computed with a new numerical method 
[T^ for likelihood ratio tests of non-nested models. If 
EWPOs and signal strengths are included in the fit we 
find p-values of 2.0 • 10~^ and 0.74 • 10~^, respectively, 
which means that the SM4 is excluded at the 3(t level. 
While this result is obtained for Dirac neutrinos, it will 
change only marginally for the case of Majorana neutri- 
nos with two fourth-generation mass eigenstates v^, v^: 
the fit to the signal sthrengths will return the same in- 
visible Higgs width, now corresponding to the sum of 
the four decay rates r{H — >• 1^4,51/4,5). A marginal dif- 
ference occurs once the EWPOs are included: choosing 
the v^-v^ mass splitting such that the eigenstate with 
the larger SU(2) doublet component becomes heavier, 
one can slightly improve the quality of the electroweak 
fit. The improvement is negligible, as indicated by the 
shallowness of the minimum of the SM4 x^ function in 
Fig. [1] While the SM4 is under severe pressure, a sequen- 
tial fourth generation may still be viable in conjunction 
with an extended Higgs sector [I31 - H7] . 
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